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Abstract 

We study a one-dimensional hamiltonian chain of masses perturbed by an energy conserving noise. 
The dynamics is such that, according to its hamiltonian part, particles move freely in cells and interact 
with their neighbors through collisions, made possible by a small overlap of size e > between near 
cells. The noise only randomly flips the velocity of the particles. If £ — > 0, and if time is rescaled by 
a factor 1/e, we show that energy evolves autonomously according to a stochastic equation, which 
hydrodynamic limit is known in some cases. In particular, if only two different energies are present, 
the limiting process coincides with the simple symmetric exclusion process. 
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1 Introduction 



Fourier's law asserts that the flux of energy in a bulk of material is proportional to the gradient of 
temperature: 

J = -k{T)VT, 

where k is the conductivity of the material. While this phenomenological law is widely verified in 
practice, its derivation from a microscopic hamiltonian dynamics still remains a very challenging 
question [2]. 

One actually knows that integrable hamiltonian systems usually violate Fourier's law, as it is the 
case for ordered harmonic chains 8 , for the Toda lattice with equal masses [9][l], or for a system of 
one-dimensional identical particles interacting only through collisions (see Remark [l] in Subsection 



2.3 1. In these three examples, ballistic transport of energy is observed. Moreover, the conductivity 



of one-dimensional chains of oscillators conserving the total momentum (no pinning) is generally 
expected to diverge [^. 

On an other hand, Fourier's law have been derived starting from some purely stochastic dynamics, 
such as exclusion processes. In these models, particles have fixed position on a lattice (pinning) and 
exchange their energy with their neighbors according to the value of some random variable. The 
simple symmetric exclusion process (SSEP) is one of them, and constitutes an easy example where 
the heat equation can be recovered from a microscopic dynamics 15:. To obtain a first hamiltonian 
derivation of Fourier's law, it seems thus desirable to start with hamiltonian systems which might 
look as close as possible to such stochastic models. In the recent years, many progress have been 
made in this direction. 

Let us consider a periodic lattice of A'^ masses, each of them being confined in a region of the 
space by a pinning potential. We assume that the strength of the interaction between near atoms is 
controlled by a small parameter e > 0. When e = 0, there is no interaction, and we suppose that the 
dynamics of each single atom has good mixing properties (hyperbolic dynamics). Then, when the 
interaction is turned on (e > 0) but is still very small, energy should fiow between atoms at a much 
slower rate than the rate at which every atom reaches its own equilibrium for a given energy. One 
expects therefore the evolution of energy to be similar to the one of a stochastic process. 

Starting from there, it has been proposed to derive Fourier's law in two steps, as described by 
Gaspard and Gilbert [4]. First, one tries to show that the energy of the atoms becomes an autonomous 
Markov process in the limit e — )■ 0, if time has been rescaled by a factor e~" for some a > 0, in order 
to obtain a non-trivial limit (but no space rescaling). The process obtained after this first step 
corresponds to a mesoscopic description of the material. Second, one hopes that the hydrodynamic 
limit of the mesoscopic process can be shown to coincide with the heat equation. 

The first part of this program has by now been carried out rigorously when the interaction is of 
the form eV, where V is a smooth potential coupling nearest neighbor atoms. Liverani and Olla [7] 
first considered the case where the uncoupled dynamics is made of a hamiltonian part perturbed by 
a stochastic noise, modeling an hyperbolic dynamics. Dolgopyat and Liverani [3] were then able to 
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obtain a similar result starting from a purely hamiltonian dynamics. One has however not yet been 
able to identify rigorously the hydrodynamic limit of the obtained mesoscopic processes. 

The aim of this paper is to investigate the case where the interaction between atoms takes place 
through elastic collision rather than a smooth potential. In this context, the parameter e represents 
the size of a small interaction zone where near particles can collide, and the limit e — > corresponds 
thus to a limit of rare interactions, rather than weak interactions. The model we consider is one- 
dimensional, and the lack of hyperbolicity of the uncoupled dynamics is supplied by a stochastic noise 



(see Subsection 2.1 1. 

The limiting process we obtain is described by ( 2.4|2.5 \ below, and belongs to a class to which 
presents probabilistic tools can be applied in order to derive the heat equation in the hydrodynamic 
regime, at least when the number of energies reachable by the system is finite 5 . In the special 
case where only two different energies are present, we even recognize the generator of the SSEP. The 
mesoscopic description of the dynamics is thus so simple in this model, that it leaves us some hope 
to go beyond the slightly artificial two steps procedure followed here. 

The rest of this paper is organized as follows. In Section [2j we describe the dynamics at a non 
technical level, we state our result, and give some notations to be used throughout the text. A 
detailed definition of the dynamics together with some elementary observations is to be found in 
Section |3] Section |4] contains the proof of Theorem [l] admitting Lemma [T] which proof is deferred 
to Section [5] The proof of this lemma constitutes actually the core of the argument. Section [5] deals 
with the convergence to equilibrium for the uncoupled dynamics. Finally, the two figures of the text 
are collected in Section [T] 



2 Model and result 
2.1 The model 

We consider a system of A'^ classical point particles of unit mass, each one evolving in a one-dimensional 
cell of size one (see figure [l] in Section [7|. Let us first describe the dynamics when all the particles 
move independently from each others. For this, let us associate a Poisson process (a "Poisson clock") 
to each of them. The A'' Poisson processes are independent but have the same parameter A > 0. A 
particle moves then freely in its cell, meaning that it travels at a constant speed and is elastically 
refiected at the boundaries. In addition, each time its Poisson clock rings, the sign of its velocity is 
reversed. So the Poisson processes serve to model a chaotic dynamics inside the cells. Observe that 
the kinetic energy of every particle is a conserved quantity under this uncoupled dynamics. 

We then introduce some interaction between the particles. Let < e < 1/2, and put the cells in 
a one-dimensional row, as depicted in figure [l] in Section [7] so that there is an overlap of length e 
between near cells. If two nearest neighbor particles hit each others, they undergo an elastic collision, 
resulting in a exchange of their momenta. This is the only interaction ; in particular, particles do not 
hit against the boundary of neighbor cells. This dynamics still preserves the total kinetic energy of 
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the system. 

To fully specify the dynamics, one finally needs to give the initial positions and velocities x = 
(gi, . . . ,qN,Pi, ■ ■ ■ ,Pn)- One supposes the initial positions to be such that particle k is at the left of 
particle k + 1 for f < fc < A'' — f . The dynamics is so that this condition holds then at all further 
times. If qk G [0, f] for 1 < fc < A*', this means qt+i > — 1 + e for 1 < fc < — I. One assumes 
also that no initial velocity is zero. For i > 0, the state of the system is then entirely given by the 
value of the Markov process 

X%x, t) = {q\{x, t),..., qlf{x, t),pl{x, t),... ,Pn{x, t)), 

where we have emphasized the dependence on e. 

2.2 The result 

Let < 61 < ■ ■ ■ < ejv' < +oo be the initial energies, where A''' is an integer smaller or equal to A^. 
These are thus numbers such that G {ei, . . . ,ejv'} for 1 < fc < A'. The dynamics is such that 
{ei, . . . ,ejv'} constitute also the set of all energies reachable by the system for all times. One wants 
to describe how energy flows between the particles in the limit e — >■ 0. Since one suspects the rate of 
collisions to be proportional to e (see Remark [s] in Subsection |2.3[ ), one needs also to rescale time by 
a factor in order to obtain a non-trivial limit. More precisely, we seek for the limit distribution 
of the time rescaled kinetic energy 

as e — >• 0. 

We need a couple of extra definitions to express our result: 

1. Let ©([0, 1], A) be the set of cad-lag functions on [0, 1] with value in A C M'*, and observe that 
actually G O([0, 1], {ei, . . . , ejv'}'^). 

2. Let 7 be a function on given by 

7(a, b) = ^ max{v^, V2b}. (2.1) 

3. Let Q"' be the physical space accessible to the particles: 

= {(gi,...,gjv) G [O.l]'^ : qk+i > qk - 1 + e for 1 < fc < AT - 1}. (2.2) 

4. A number ^ G R is called diophantine if there exist constants C,/3 < -|-oo such that, for every 
p/g G Q, one has 

With the definitions and notations introduced up to now, one has 
Theorem 1. Let ei(0), . . . , ejv(O) G {ei, . . . , e^v'}. Let us assume that 



1. for 'i < j k < N' , the numbers ^/ejjek are diophantine, 
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2. the law of X'^( • ,0) is the only probability measure which density is proportional to 



XQ^ill, • ■ • ,giv) •X(ei{0),...,e„{0))(Pl/2, . . . ,Piv/2). 



Then 



1. as e ^ 0, the distribution of the process tends weakly to the distribution of a process 



£ = {£^,...,£n) G O([0,l],{ei,..., 



ejv'} j, 



2. £ is the unique Markov process which law P solves the equations 



JV-l 



atP(ei,. . .,eN,t)= ^ 7(efe,efe+i)fp(. . . , efc+i, Cfc, . . . , t) - P(. . . , efe, efc+i, . . . , f)j , (2.4) 



Remark 1. The dynamics, with or without noise, should actually becomes completely elementary if 
the particles were not confined into cell^ Indeed, in this case, the system can be seen as a collection 
of N independent particles, since they just exchange their momentum when colliding. 
Remark 2. The condition on the initial measure is much more restrictive than needed. A close 
look at the proof shows that, with some extra work, one should be able to start from a measure with 
bounded density with respect to the positions. Our choice mainly allows us to avoid an extra initial 
step that should be needed to apply Lemma [l] below. 

Remark 3. It is not so obvious that time has to be rescaled by a factor in order to obtain a 
non-trivial result (the correct scaling factor is in [3] and [7]). So let us try to motivate this a 
little bit, following the heuristic of |4]. It is enough to concentrate on the case of only two particles 
(A^ = 2), with two energies ei < 82. When e becomes very small, the rate of collisions between these 
particles should tend to zero as well. So, since the dynamics of uncoupled particles have good mixing 
properties, the distribution of position and velocity of a given particle with energy e G {si, 62} should 
be well approximated by the equilibrium measure of an uncoupled particle of the same energy. It is 
easily checked that this measure is given by 



Assuming that, at a given time, the energy of particle 1 is ei, and the energy of particle 2 is 62, it 
seems then reasonable to identify the instantaneous energy exchange rate as 



fe=i 



P(ei, . . . ,ejv,0) = X(ei(o),...ejv(o))(ei, . • . ,e]v). 



(2.5) 



2.3 Remarks 




■7 = lim — P f particle 1 hits particle 2 in a time interval n ) 
j7-t-o r) ^ ' 




^ we thank J.L. Lebowitz for having recalled us this fact 
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where one has taken into account the fact that particle 1 should be at the left of particle 2. One 
computes that 

7 = — ^ — = -max{V2ei, V2e2}, 



so that indeed 7 will coincide with 7 as defined in (2.1 1 if time is rescaled by a factor e~^. 
Remark 4. The presence of a diophantine condition on the velocities may appear as a surprise in 
the present context. To see where it comes from, let us again take the case N = N' — 2, and let us 
assume that one velocity has value 1, and the other 2, violating thus the diophantine condition. The 
dynamics of a single particle is such that, for any fixed time interval, in the microscopic time scale, 
there is a strictly positive probability, independent of e, that the Poisson clock of the particle does 
not ring during this time, so that its trajectory is in fact deterministic. Now, let us suppose that the 
two particles collide, and let us take the "typical" case where the collision does not occur too close 
to the borders of the overlap region between the two near cells. It can then be seen that, no matter 
how small e is taken, they will collide once more with each other, in a time interval of order 1, if they 
travel along the deterministic trajectory. This happens thus with a strictly positive probability. 

However, the basic assumption behind our theorem is that successive shocks occur "randomly", 
and are typically spaced by time intervals of order , in the microscopic time scale. This needs 



to be so for the coefficient 7 to be given by (2.1 1 in (2.41. But we see that, if the quotients of 



velocities are rational like here, successive shocks are actually quite correlated, and tend to occur 
in clusters. It seems clear however, that these clusters should themselves occur at random and be 
typically spaced by time intervals of order e^^, so that we conjecture Theorem[l]to still remain valid 
without diophantine condition, but with a smaller coeffeicient 7 since recollisions tend to cancel the 
transfer of energy. We have not been able to show this. 

2.4 Notations 

Constants. The number called constants in this text may depend on the following parameters of 
our problem : the number A'^ of particles, the number A*'' of different energies, and the value of the 
smallest and the largest energy, namely ei and ejv' . They never depend nor on e, nor on time. 
Time discretization (rj). We will find convenient to divide time in small intervals. Throughout 
the text, we will denote these time steps by 77. One always assume rj < e. Its exact value matters 
few, and one can always think that rj « e. 

Deterministic dynamics (~). Objects with a tilde refer to the deterministic dynamics. So e.g. 
X'^ is the trajectory of our process driven by the deterministic dynamics. 

The phase space. Points of the phase space $ are written x — (q, p) — (qi, . . . , qN,Pi, ■ ■ ■ ,Pn)- 
Moreover, one writes e = (ei, . . . , Cn) ~ {p\/2, . . . ,p%/2). 



Norms. Let G be the Gibbs measure defined in Subsection 3.3 below. The norm of the space L''(<I>, G) 
is written || ■ ||p (1 < p < +cxd). If is a measure on one writes ||/i||i = sup^,.||^,||^<j^ \ij.(u)\. 
The space 5*^. Functions on $ of the form u(q, p) = XQ' (q) • ^'(e) are said to belong to the space 
. A measure p on $ which density w.r.t. G belongs to , is itself said to belong to 5". 
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3 The dynamics 



3.1 Definition 

Let e > 0, let TV > TV' > 1 be two integers, and let < ei < ■ • ■ < Sjv' < +oo be given energies. Let 

P = {-^/2ei:,V2er,...,-V2^,^2^} 

be the set of velocities, let 

$ = [0,1]^ xP^ 

be the phase space, and let 

= Q'x P^, 



with Q" defined by 1 2.2 1, be the subset of <I> where the dynamics actually takes place. Let x = 
(q. P) = {gi, • • • I qN,Pi,- ■ ■ ,Pn) e "l" and, for t>0, let 

X^{x,t) = {ci{x,t),p^x,t)) = {ql{x,t), . . . ,qlf{x,t),pl{x,t), . . . ,pN{x,t)) 

be the stochastic process on "I> such that X'^{x,0) = x and that, dropping out the superscript e for 
clarity, 

dqk =Pkdt, (3.1) 
dpk = - 2pfc(dNfc + dN'fc + dN'fc) + (pfc-i - pk) dNl_i - {pk - pk+i) dNl (3.2) 

for 1 < fc < TV, where 

1. the processes are independent Poisson processes with identical parameter A > 0, 

2. the processes N\ are cad-lag processes such that N5j.(t+) — N|.(t_) G {0, 1} for t > 0, with 

Ni(t+) -Ni(i_) = 1 if and only if gfe(t_) = and pfc(t_) < 0, 

3. the processes are cad-lag processes such that N's.(t+) — N'j.(t_) G {0, 1} for t > 0, with 

Nfe(i+) - Nfe(i_) = 1 if and only if qk{t-) = 1 and Pk{t-) > 0, 

4. the processes are cad-lag processes such that N^(t+) — N^(t_) £ {0, 1} for t > 0, with 

N^(f+) -N"fc(t_) = 1 ifandonlyif gfc+i(t_) = gfe(t_) - 1 + e a.nd pk+i{t-) < pk{t-) 

for 1 < fc < TV - 1, and = 0. 

Here, f{t+) and f{t-) represent respectively the left and right limit of a cad-lag function / at time t. 

Equations ( 3.1|3.2 \ describe the evolution of a particle in a TV-dimensional billiard, which projec- 
tion on the (gfc, gfe+i)— plane is depicted on figure [5] in Section [t] It is checked that ( 3.1|3.2 \ admits 
a unique solution given almost any initial condition w.r.t. the measure G" defined in Subsection 



3.3 



below. Indeed, looking at figure |2] the deterministic solution (when none of the Poisson clock rings) 
can be constructed at hand for a fixed finite time interval. The general case can then be obtained, 
noting that the Poisson processes have almost surely finitely many jumps in any time interval. 
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Observe that the condition pk+i{t-) < Pkit-) in the definition of NJ: is such that particles may 
move from <I> — $^ to but may not escape from <!>''. The definition of the dynamics outside is 
absolutely irrelevant for the statement of Theorem [T] but our choice turns out to be convenient for 
its proof. 

3.2 Evolution semi-group and generator 

Let e > 0. For t>Q, the operator "P^'* acts on bounded functions u on $ according to the formula 

V^'\{x) = E{uo X'{x,t)). 
Its adjoint p*^*'* acts on measures on and is defined by means of the relation 

Let now u be some smooth function on $, and let v = v{x,t) := ■'P^u{x). The evolution 

equations imply that ( 3.1|3.2 \ that v satisfy the boundary conditions 



v{. . . ,qj, . . . ,pj, . . . ,t) = v{. . . ,gj, . . . , -pj, ...,t) if = oi = 0, (3.3) 
'"{■ ■ ■,<lj,<lj+i, ■ ■ ■ ,Pj,Pj+i, ■■■,t) =v{.. .,qj,qj+i, . . . ,pj+i,pj,. . .,t) if qj+i = - 1 + e, (3.4) 

the first of these relations being satisfied for 1 < j < N , and the second for 1 < j < N — 1. Moreover, 
it solves the following differential equation in the interior of $: 

JV JV 

dtv = ^pkdg^v + X^{v{...,-pk,...)-v{...,pk, ■■■))■ (3.5) 

k = l k=l 

3.3 Invariant measure 

The probability measure is defined as the uniform probability measure on 



pSP' 



where 



IP^^I = (2iV')'^ and IQ'I = / dq. 
When e = 0, we will just write G instead of G''. One checks using (3.51 together with the boundary 



conditions ( 3.3|3.4 1 that the measure G^ is invariant under the dynamics, meaning that 7-"^*'*G'^ = G"^ 



for every t > 0. Since energy is conserved, this is not the only invariant measure. In fact G"^ 
corresponds to a Gibbs measure at infinite temperature. 

4 Proof of Theorem [I] 

Let < eo < 1/2. One first shows that the family {£'^)^<^g is tight in T>{[0, 1]). One then proves that 
the law of £'(t) converges weakly for every t G [0, 1] to the law P{- ,t) characterized by ( 2.4|2.5 1. One 
finally establishes that the limit process £ is markovian. 
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4.1 Tightness 

We prove here that the family )e<eo is tight in ©([O, l],{ei, . . . ,ejv'}^), meaning that for every 
(5 > 0, there exists a compact set Kg C 25([0, 1], {ei, . . . , e]v'}^) such that P''{Ks) >l — 5 for every e G 
]0, eo], where P*^ is the law of S'^. For n > 1, let K„ be the set of functions in ©([0, 1], {ei, . . . , e]v'}^) 
having at most n jumps. Since {si, . . . , Sjv'}'^ is finite, the sets Kn are compacts. Let now 5 £]0, 1[, 
and let us show that there exists n{5) > 1 such that P^{K„(^s)) > 1 — 5 for every e g]0, eo]. 

The energy of the particles can only change due to a collision with a neighbor, and so it suffices to 
establish that there exists a sequence {c„)„>i, with c„ — >■ as n — >■ oo, such that, for every e g]0, eo] 
and 1< fc < iV- 1, 



P. 

where ^ G 5*^ is the law of X'^{ • , 0). 



dNlis) > < c„, (4.1) 



One would like to replace the integrand in this expression by an explicit function on the Markov 
chain X^. If the dynamics were deterministic (no Poisson processes), one should have 

r dNUs) < C r ■XY,i.,^)oXtds (4.2) 

Jo Jo 

with r; > as small as we want, and 

Yfe(e,77) ;= {a; G $ : jgfc+i - + 1 - e| < r?}. 
But, for e > one finds 77 = i]{e) small enough such that the event 

A(e,77) — {oj: max (Nfe(i + 77) - Nfc(t)) + (Nfc+i(i + r?) - Nfc+i(t)) < l} 

te[o,<!-i] 

has a probability as close to 1 as one wishes. 



Therefore, taking rj small enough for a given e, one may replace (4.1 1 by 

P^(^' dNUs) > n|A(e,r,)) < c„. (4.3) 



Doing so, inequality ( |4.2[ | may now be applied (with a larger constant) even in the presence of the 
Poisson processes, and, by Markov's inequality, one gets 

= Cn~\^ dsy^P''*(77"\xYfe{e,^))dp 
< Cn"^e"\G'(7)"^XYfc(£,,,)) < Cn"\ 

where one has used the fact that < C to get the second inequality, and the fact that ^{r]'^ . XY^(e,r])) 
0{e) for any value of ri g]0, e] to get the last one. □ 

4.2 Convergence at a given time 

We show here that the law of£'^{t) converges weakly for every t G [0, 1] to the law P( • , t), characterized 
by ( 2.4|2.5 1. Observe that, for given e > 0, there is an obvious identification between functions on 



{ei, . . . ejv'} and measures in 5'^, that will be used in the sequel. 
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Let P'^{ - ,t) be the law of X'^{t) (so not the law of £"{1) as it was in the previous subsection). We 
actually will prove that, for any t G [0, 1], there exists a probability measure P{ - ,t) G 5", such that 
P'^( ■ , t) converges weakly to P( • , t) as e — > 0, and then that P{ - ,t) solves 
our result. Indeed, the distribution of £'^{t) may be written as J^^ P''(d2;, t|e), and so, as e — > 0, the 
distribution of f^(t) also will converge to P{ - ,t). We will use the following lemma, which proof is 
deferred to the next section. 

Lemma 1. Let v be a function on {ei, . . . ,ejv'}^ corresponding to a probability measure m S" for 
every e £]0, eo]- For r > small enough, and for < e < r^, 

N-1 JV-1 

P*'''-''^P = (l-r ^ 7(efc,efe+i)) P + r ^ 7(6^, e^+i) P(. . . , efc+i, e^, . . . ) + 0(r')M.,., 

1 k = l 

where P £ is the measure determined by v, and where fx^^T is a measure on $ such that ||/^E,r||i < 1. 

Let now t G [0, 1]. Let e g]0, eo[, and let be the largest integer such that r :— t/n^ > e^^^. One 
has P'{-,t) =-p"''"'*P'(-,0). Lemmajl] applies if e is small enough, and one has 

p"''"'*p(-,o) = r"''''"'^p{- ,0) = p"'^"'<"^-''"(Pi + o(nr') Mi...r), 

where Pi G 5"^ is the measure 

N-1 N-1 

Pi ~ (^1 - T ^ 7(efc,efc+i)j P'(- ,0) + r ^ 7(efc,efc+i) P''(. . . ,efc+i,efc,. . . ,0), 
fc=i fc=i 

and ^i,e,T a measure such that ||/ii,e,T||i < 1. Observe that Pi is a probability measure if r is small 

enough, what we assume. Since ||7''^*'°At|| 1 < for every measure fj, and every s > 0, one may 

iterate this and get 

Tig 

P"'^"''P^(-,0) = P„. +C'(n7')2Afc,,,,, 

fc=i 

with Pn, a probability measure in S'^, and fik,t,T measures such that ||/ifc,E,T||i < 1. 

Let us first show the existence of a probability measure P( ■ , t) G iS'' in the closure of {Pn,)t<tQ for 
the II ■ ||i-norm. For this, let Vn, be the function on {ei, . . . , ^n'}'^ determining the measure P^^ G 5", 
and let then P„^ be the unique measure in S'^ determined by Vn^- Since \\P„^ — PnJIi < Ce^, a point 
in the closure of (Pnj)e<Eo is also in the closure of (Pnj)e<eo and is a probability measure. But, 
since 5" is a finite-dimensional space, and since (P„,)e<£o is a bounded set for the || • ||i-norm, the 
announced measure P{ - ,t) indeed exists. 

Let us next show that P( ■ ,t) solves 

P{-,t) = limP'=*''"'*P(-,0) (4.4) 

e— )-0 

for the II - ||i-norm, and so 

^.^ P(.,t + r)-P(.,t) ^ ii^iiim(P"--^--Id)P"--^*P(.,0) 

T-^0 T T->0 r e-»0 ^ ' 

= lim if lim(P"''"'"-Id)P( ■,*)') 

+ l™o 7 ( 1^0 (^"''"^ - ^'^) r *""*P( • , 0) - P( ■ , t))) . 



2.4|2.5l. This will imply 



2.4|2.5 I . Up to a subsequence, one can write 
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The second term in the right hand side of this equation is zero since 

||(73">^-'- _Id)(p"-'="'*P(.,0)-P(-,t))||^ < 2||p"''="''P(.,0)- P(-,t)||^ ^0 
as e — >■ 0. By Lemma [l] instead, the first term gives 

lira - ( lim (p"-^"'- - Id) P( ■ , t)) = X] 7(efc, efc+i) (P(. . . , efc+i, e^, . . . ,t) ~P{. . . , e^, e^+i, . . . , t)) , 

fc=i 

which estabhshes the result. □ 

4.3 S is a Markov process 

We show that, for any = to < ii < • ■ • < in < 1 and any eo , ei , . . . , e„ G {si , . . . , ejv } , for n > I , 

P(£t„ = e„jft„_j = e„_i, . . . ,fo = eo) — P(£t„-t„_i = e„ |£o = e,i_i) . 

This is immediate for n — 1. We show the case n = 2, the generalization being rather straightforward. 
For e > 0, since is a Markov process, one may write 

P{£t2 = (^2\£ti = ei,£:o = eo) = P^^ {£t2-ti = 62) 

where /x'^ is an initial measure on <I> given by 

fJ." = V"'^ *Vo(-le = ei) with /iq = G^(-|e = eo). 

Our claim will now be established if we show that H/i*^ — G'(-|e = ei)||i — )-Oase— >-0. But this follows 



from the facts that WV"-'- *Vo-P( • , ti)||i as e by ( [XI] ), that P( • , ti|e = ei) = G( ■ |e = ei) 
since P(- G S°, and that ||G - G'||i = C(e^). □ 



5 Proof of Lemma [l] 

The proof of Lemma [l] is conceptually very easy. We start by sketching the method in Subsection 
|5.1| In Subsection |5.2[ some crucial elementary facts are collected. We finally prove Lemma [l] in 
Subsection 15.31 

5.1 Method 

Since P G 5^, this measure is invariant under the uncoupled dynamics, up to an error of order for 
the II • ||i-norm. We recall that e < by hypothesis. Writing the result of Lemmajljas 

one sees that one needs to evaluate the difference between the evolution of the coupled and the 
uncoupled dynamics in the first order in r. The following identity, known as Duhamel's formula, is 
suited for this: if 5* and T are two bounded operators on some space, one has 

n 

T"-s" = ^^*=~'(^-s)S'""^ (5.2) 
fc=i 



11 



Iterating this formula to replace T by 5* , one finds 

n n k — 1 



fc = l 3 = 1 



(5.3) 



k=l 



Let now rj be some small time step (see Subsection 2.4 1, let S = 'P"*-'', let T = j}^*'^^ and let 



n be such that e = n?). At this point, it could have seemed more natural to work directly in a 



contirmous-time setting. In that case, the sums appearing in (5.2 1 and (5.31 becomes time integral. 



whereas T — S should become the part of the generator corresponding to a collision between particles. 
The problem is that this term is hard to make explicit, since it is only present through the boundary 
conditions ( 3.3|3.4 l, and that is why we decided to divide time in arbitrary small intervals. 



All what needs to be done now is to show that the first term in the right hand side of (5.3 1 



coincides with the first term in the right hand side of (5.1 1 up to an error of order e, and that the 
second term is O(r^). The identification of the first terms is not hard since one has a good control on 
the uncoupled dynamics (see Section l6||, i.e. on S = 75°*'", and since one can estimate the difference 
T — S = p"*'"^ — •po*'i for ri small enough, as explained in Subsection 



5.2 



below. 



The second term represents recoUisions. One needs to see that, in a macroscopic time of order r, 
the probability of two collisions or more for the same particle is C'(t^), as it should be if they were 
happening at random. This turns out to be a bit heavy to show, mainly due to the persistence of a 
deterministic trajectory for microscopic time intervals. The diophantine condition on the velocities 
has here to be used. 



5.2 Crucial elementary estimates 

Let e > and let e > > be a small time interval. Let us define some subsets of For this, let 

p ;= 2max{V2ei, . . . , \/2ejv'}. 



Then, for 1< fc < iV - 1, 



Zfe = {a; G $ : Igfc+i - gt- + 1 - e] < prj}, 
CoUfc ={x €<!>■. Nl{x,r,)~Nl{x,0) = 1}. 



3.1 



(5.4) 
(5.5) 

for a deterministic 



Here N^(3;,s) is the value at time s of the process defined in Subsection 
trajectory starting at x. We recall that the use of a tilde always refers to the deterministic dynamics 
(see Subsection |2.4[ ). One then sets 

JV-1 JV-1 

Z = (J Zfc, Coll = y CoUfc. 

fc=i fc=i 

First, for 1 < fc < A'^ — 1, the characteristic function of the set Zk depends only on the variables 
Qk and Qk+i- A computation yields 



Cerj < / XZfc dqkdqk+i < / Xz dq < C' erj, 

i[0,l]2 Jq 



(5.6) 
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and, for 1 < j / fc < iV - 1, 



/ 



Xz, .Xz.dq = 0{erif. (5.7) 



/q 

Second, for 1 < fc < A'^ — 1, let us define 

Wfc = G $ : 1 - (7fc < pry and \qk+i - + 1 - e| < pr;}, 
Wfe = {x G $ : gfc+i < P77 and \qk+i - + 1 - e| < pr?}. 

The deterministic dynamics is such that 

Coll C {x G <!> : pfc > pk+i and < qt+i - gfe + 1 - e < (pfc - Pfc+i)??} U Wfc U Wfc, (5.8) 
Coll D {a; G $ : Pft > Pk+i and < qt+i - + 1 - e < (pfc - Pfc+i)r?} - (Wfc U Wfc). (5.9) 

So CoUfc C Zfc and thus Coll C Z, and the characteristic function of CoUfc only depends on the 
variables Xk and Xk+i- A computation based on ( 5.8|5.9 1 gives 



/ XCoiifc dgfcdgfc+i = XK+(Pfc -Pfc+i) • (Pfc -Pfc+i)'7e + ©(r?^), (5.10) 



for any values of the impulsions. It follows moreover from (5.71 that, for ^<j^k<N — 1 



XCoiij • XCoUfc dq = 0{eri) (5.11) 

for any values of the impulsions. 

Let us finally estimate the difference between the evolution of the coupled and uncoupled dynamics 
over a time interval rj. Let us define the operator 

L ■- 1(7'^'" -pO'"). (5.12) 



One has 



and 

Lemma 2. There exists a constant C < +oo such that, for all rj > small enough, and for all 
uG L°°($), 

Lu = Lu + Xz n 
where u' G L°°(<1?) is such that ||u'||oo < C||u||oo. 

Proof. Let a; G If a: ^ Z, then X\x,ri) = X°(a;,r?) for every realization of the Poisson processes. 
Therefore 

Lu — Xz . Lu. 

Next, the probability that one of the Poisson processes have a jump in a time interval of length 77 is 
itself 0{Tf) so that 

V'^'uix) = E{u o X\x, r,)) = uo X'{x, r;) + O{rj\\u\\oo) = 'P''''u{x) + 0{r]\\u\\^). 

This last formula still holds also in particular for e = 0, and therefore Lu = ^^Lu + Xzit' with 
||m'||oo < C \\u\\oo for some C < +00. □ 
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5.3 Proof of Lemma [T] 

1. Preliminary simplifications. We will give the proof in the case TV' = A'^, indicating at the 
very end the small needed adaptations to bring for dealing with the case TV' < TV. The space S'^ is 
generated by functions of the form 

Xefc^,...,efc„(ei,...,ejv).XQ-(q), 

where fci, . . . , fcjv is a permutation of 1, . . . , TV. By linearity, it suffices to consider the case where 
P is proportional to a function of that form. Moreover, to simplify notations, we will assume that 
{ki, . . . , fcjv) = (1, . . . , TV). Let then 

P' := ^-Xei e„(ei, . . . ,e]v) .XQ(q) = ^ • Xei,...,e„ (ei, . . . , eiv). (5.14) 

This measure is close to P when e — )■ 0: 

||P-P'|li = O(e'). 

Because HP'*'*/^!!! < for every /i and every t > 0, one may therefore show the result with P' 

instead of P. This will be advantageous since P' is invariant under the uncoupled dynamics. 

2. Applying Duhamel's formula. Let r > 0, let e G]0,r^[, and let 77 e]e*,2e'*[ be such that, 
for some integer n > 1, one has 

e~^T = rjn. (5.15) 
Applying Duhamel's formula (15. 21) and noting that p^'-'p' — p' for every t > 0, one gets 



with L the operator defined in (5.121. One would like to replace L* by L* , with L defined in (5.131 
One has 



|P"'<'=-''''(L* -Z*)P'||i < ||(i* -Z*)P'||i = sup P'{{L-L)u) 

":|klloo<l 



< sup max / j (i — L)it(q, p)| dq = 0{eri). 

u:\\u\\^<l PSP" Jq 

Here the inequality follows from the fact that P' is a probability measure independent of q, whereas 
the last equality is a consequence of Lemma [2] and ( |5.6[ ) . Therefore 

n 
fc = l 

where ^0 is a measure such that 

IImoIIi = 0(??n.e77) = 0(r'). (5.17) 
One then uses Duhamel's formula once more to replace ■p'^*'('=~i)'' by ■p''*'('=~i)'' in (5.161: 



fe=l fe=l j=l 

:= P' + Ml + M2 + Mo- 
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It is not clear at this point that we can replace the remaining operator L* by L* up to a negligible 
error, and so we directly proceed now to the analysis of /xi and ^2- In a first step we will show that 



Hi is given by (5.231 below, whereas in a second step, one will prove that ||/i2||i = 0{t ). Thanks to 



the bound (5.171 on H/ioHi, this will finish the proof. 



3. Analyzing /ii. We here show that /^i can be expressed as in (5.231 below. Given a measure 
/i on $, we write /i(l|e) the measure in iS*^ explicitly given by 



(5.18) 



where the factor 1/2^ appears since we want to see /i(l|e) as a measure on <!>, and not as a measure 
on {ei, . . . , Biv}^. Because P'^*'*/i(lje) — ^(l|e) for any /i on "1> and t > 0, and since rjn = e^^r, one 
has 

n 

fii = e-V.L*P'(lle) + 77^P°*'"=-')''(L*P' - Z*P'(l|e)). 



(5.19) 



Let us first compute L*P'(lje). It follows from ( |5.18 1 that 
Z*P'(l|e 



Now, from (5.131 and then (5.111 



^'(LXe) = i (P'(XCon . ^^'"Xe) - P'(XColl . ^'''''Xe)) 

1 

= ^ E (p'(xcoii, .P^'^xe) - P'(xcoii, .P"'"xe)) + o{n). 



On the one hand, from the definition (|5.5|), one has 



with 



^"''Xe(a:) = Xe(^'(2;,77)) = X^.e(x) for 2: G CoUfc - (J CoU, , 



(Tfce = (ei, . . . ,efe+i,efc, . . . ,ejv). 
On the other hand, P^^'^Xe = Xe since the uncoupled dynamics preserves the energy of individual 



particles. Therefore, using (5.111 once more, 

JV-1 



= ^ E (P'(^Coll,.X<..e)-P'(xColl,.Xe)) + 0{V) 



From the particular form ( |5.14[ | of P' and from (5.101, one obtains 

P'(XCoiifc .X<Tfce) = E X{ei,...,e„)(p') •X<Tfee(p') ■ / XCoiifc (q) dq' 



p'6P' 



IE 

2JV 



E X{ei,...,e„)(p') •X<Tfce(p') •XK+(p'fc -P'fe + l) • (Pfc -P'fc + l) + <^('7^) 



p'6P' 



X(ei,...,e„)(o-fce) ^ XR+ip'k - p'k+i) ■ {p'k - p'k+i) + Oirj"^) 



»7e-X(ei,...,ejv)(o-fce) .7(efc,efc+i) + 0(77 ), 
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where 7 is the function defined by (2.f I. A similar computation shows that 



P'(XCoii;, -Xe) = r;e.X(ei,...,e„)(e) .7(efe,efe+i) + ©(r;^), 
so that one concludes that 

N-l 

L*P'(ile) ^ 7(efc,efc+i) . (xei,...,e„ (c^fce) - Xei,...,e„ (e)) +£>(??) 



e^7(efc,efc+i) • (P'K-e) - P'(e)) + 



(5.20) 



We then need to bound the second term in the right hand side of (5.191. Let us first observe 



that, by the particular form ( |5.14[ ) of P', by the definition ( |5.13[ ) of L, by the fact that Coll C Z, and 



finally by ( 5.6 1 



L*P'||i < sup max / jl/?t(q, p)|dq < max — / xz(q)dq = C(e). (5.21) 

u:||u||^<l peP« Jq pGP« V Jq 



The uncoupled dynamics is studied in more details in the next section ; it follows from ( 5.21 \ and 



(6.81 there that 



V 



n 71 



(5.22) 



where c denotes some strictly positive constant. Putting (5.20 1 and (5.221 in (5.191, one ends up with 

iV 



1 = T^7(efc,efc+i) . (P'(crfce) - P'(e)) + fi'i, 



(5.23) 



with ll/^'illi = 0(t'^), since 77 < e < r^. 



4. Bounding fi2- One here establishes that ||/i2||i = 0{t'^). Since < ||/i||i for every 

s > and every measure fi, one has 



1 < (e"V)77^ ||L*P''*'(^-i)''L*P'||i, 



(5.24) 



and, because of the particular form ( |5.14[ ) of P', 



||^.pu.,u-i)^^*p'||^ < sup max / lLP"'^^-'"'L?i(q, p)| dq. (5.25) 

So let us take u G L°°($) such that \\u\\ac < 1, and p G such that p'i = 2efe for 1 < < iV. To 
lighten some further notations, let us also define the time 



By Lemma [2j 



i = (j - 1)'7- 



(5.26) 
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with llu'lloo ~ To avoid to deal with too many constants, one will assume that actually 



|it'||oo < 1- Then, from the definition (5.131 of L, 



and, similarly. 



Therefore 



< (^XCoii • ■P'''''P''''xcoii + xcoii • ■p'^''''P°'*xcoii^ 



\Lr"'\zu\ < ^(xcoii.P'^''P"' xz + XCoii.P'''''P°'*Xz). 



|LP"''LM(q,p)|dq < - 

Jq 

2 



XCoii(q, p) . ■p'''''P°'*XCoii(q, p) dq 

+ ^ / XCoii(q,p).P°''''P''''xCoii(q,p)dq 

V Jq 

+ - XCoii(q,p) •P''''P°'*Xz(q,p)dq 

V Jq 

+ - f XColl(q,p).P'''''P°'*Xz(q,p)dq. 

V Jq 



(5.27) 
(5.28) 
(5.29) 
(5.30) 



The four terms appearing in the right hand side of this inequality are bounded in a very similar 
way, and we will only deal in detail with the first one, indicating at the end the minor needed 
adaptations to bound the three others. We write 



XCoii -P ' P ' XCoiidq < 

Q l<fe,!<JV-l 



E 



XCoiifc .P'''''P°'*XCoii, dq. 



(5.31) 



It follows from ( 6.6|6.7 1 and then (6.3 1 in the next section, that 



■p°'*XCoii, (a;i,a;i+i) = 



+ e" 



+ e" 



* • XCoii, {Xi{xi,t),Xf+i{xi+i,t)) 
J gixi+i,x'i+i,t) .xcoUi{Xi{xi,t),x'i+i)dx'i 

g{xi ,x[,t). xcoii, {x'l , Xf+i {xi+i, t)) Ax'i 



i+i 



+ / g{xi,Xi,t) .g{xi+i,x'i+-^,t) .xcoiii(x')Ax\Axij^-^, 



where we have used the notation 

uAXm 



E 

Vm=±\Pry 



udq„ 



l<m< N. 



[0,1] 



The first term in the right hand side of this expression cannot be simplified further. As far as the 
second is concerned, by the bound Coll; C Z; and the definition (15. 41) of Zj, one obtains 



XCoUi{Xi{xi,t),x'i+i) < Xli-^,i]{^{xi,t)) .X[o,pr,](|?;+i - qi{xi,t) + 1 - 
Now, g is uniformly bounded according to Lemma [S] in the next section, and so 

g{xi+i,x'i+i,t) .xcoU,{Xi{xi,t),x',+i)dx',+i = 0{f]) ■ XIi-^,i]{Qi ixi,t)) . 
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By a similar computation for the third term, and using (5.61 for the fourth one, one concludes that 



+ Oir,).xli-,,i]{q?{xi,t)) + C>(r,).X[o,.](glVi(a;i+i,t)) + Oir,e) 

4 

: ^A,;,,j(a;,,a;,+i) (5.32) 



where the index j comes from the fact that t = {j — l)r; by (5.261. 
Defining then 

B,,j,fe,;(pfc,pfe+i) ^ \ j xcoiifc .■P''''Aij,i dq, 



and going backwards in the expressions one concludes that it will be enough to show 



that, for i — 1,2, 3, 4, and for 1 < fc, / < TV — 1, one has 

n 

r;X]B,,,,fc,,(pfe,pfc+i) = 0{re), (5.33) 
which in particular will be the case if one establishes that 

Five cases will be analyzed separately: {i = 1,1 — k), (j = 1, Z = fe ± 1), (i = 1, |/ — A:| > 1), (i = 2, 3), 
{i — 4:). The diophantine condition on the velocities only needs to be used in the case {i — 1,1 = k). 



1. Bounding Bij,fe.fc. One has 



Bij,fc,fc < — — xcoiijx) .xcoiifc(Xfe(Xfc(a;,77),f),Xfc+i(Xfe+i(x,ry),t)) dq. (5.34) 
V Jq 

Let us take some x £ CoUfe. The first point to observe is that, for the deterministic realization X' of 
the process X'^, collisions involving a given particle are spaced by time intervals of order 1 at least. 
There exists therefore a constant c > such that, for t = {j — 1)77 < c, 

{xUXlix,7^),t),x2+^(Xl+i{x,v),t)) i Collfe, 
so that Bi_j_fe,fe = in that case. Let then suppose t>c, and let us just bound the second factor in 



the integrand of (5.341 by 

XCol\k{Xk,Xk + l) < X[l-2£,l]x[0,2el ((?*:, gfc + l), 



valid since CoUfc C Zfe, and by the definition ( |5.4[ ) of Z^. Let then t* > c be the smallest time such 
that 

{ql{XU^,V),n,ql+i{Xl+i{x,7^),n) e [l-2e, 1] X [0,2e]. 
There has to be two integers r,r' ^ such that 

2r- 2r' 
t* = — +0(e) = — +0(e), 
Pk Pk+i 



"I - 0{e) and = 0{e) 



and so 

Pfe+i r' I " ' I pfc r 
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As said after (5.251, p| — 2efe and pl+i = 2ek+i, and, by hypothesis, the quotients pk/pk+i and 
Pk+i/pk are thus diophantine by hypothesis. This condition guaranties the existence of a constant 
c > such that r, r' > ce~^^^, and so one has also t* > c'e"^''^ for some c' > 0. Since Bij^k,k = 
as long as t = (j — l)r; < t* , and since rj > e*, one concludes that 

Bi,,,fc,fc < = O(e'). (5.35) 

2. Bounding »7X]j=i Bi,j,fc,fc±i. The two cases, I = k + 1 and / = fe — 1, are analogous, and one will 
treat the first one only. One has 

Bij,fe,fc+i < ^ [ XCoiiJa;) .xcoiifc+i (^fc+i(^fc+i(3;,»7),t),X°+2(^fc+2(a;, dq. 

Integrating over qk and Qk+i, one gets 



Bi,i,fc,fc+i < — sup / xzfc+i (Xfc+i(Xfc+i(a::,77),f),Xfc+2(^fc+2(a;,r?),t)) d(jfe+2d(jfe+3 

^7 gfcS[l-2e,ll J 
9fc+l6[0,26l 

One may obviously insert the factor X^h^i + (1 ~ Xz^+g) ^^is integral, and split it accordingly. In 
the term containing the factor xZfc_|_2' ^'^'^ may just use the bound xz^+i (...)< 1 so that this term 
is 0(776), whereas, in the term containing the factor (1 — X^k^i)^ '^^'^ may replace by Xfe_|_2. At 

the end, 

Bi,i,fc,fc+i < Ce^ H sup / XZfc+i(Xfc+i(^fc+i(a;,»?),t),^fe+2(a^,i + J7))dqfc+2. 

"n 9fce[i-2E,i] J 
<jfc+ie[o,2£] 

For the integrand to not be 0, the two following conditions need to be fulfilled: 

gfc+i(Xfc+i(a;fe,a;fc+i,7?),t) G [l-2e,l], 
|<?fc+2(a;fc+2,i + f?) -^fc+i(i7fc+i(a;fe,a;fc+i,'7).i)| = 0{rj). 
The uncoupled dynamics is such that, for any a G [0, 1], any s > 0, 

Leh{\qt+^{xk+2,s)-a\=Oiv)) = 0(77). 

So, integrating over qk+2, 

Bij,fc,fc+i < Ce^ + Ce sup Xli-2,,i]{qk+i{Xl+i{xk, Xk+i,ri), {rj - , 
9fce[i-2E,i] 

9fc + l6[0,2El 



where one has used the definition (5.261 of t. Let us thus fix some qk (z [1 — 2e, 1] and some 
Qk+i G [0, 2e]. Summing 'Bi,j,k,k+i — Ce^ over j, one sees that only a proportion e of the terms are 
non zero, and so 

n 

rjY,B^,j,k,k+i = 0{er). (5.36) 

3. Bounding Bij,fej, |i — fcj > 1. One has 

Gi,j,k,i < XCoiifc(a::) .XCoii, (Xf(Xf(a;,?7),f),X;+i(Xf+i(a;,77),t)) dq. 



Jq 



19 



Integrating over the variables qk,qk+i, one gets, 



Bi,j,fe,; < sup 

ri gfe6[l-2e,ll 
qfc + ie[0,26] 



where dqfe,fc+i = dgi 



5rfe_idgfe+2 • • • dgjv. One may obviously insert the factor 

Xz,_iuz,uz,+i(a::) + (l - XZi_iUZ,uz, + i (a::)) 



into this last integral, and split it accordingly. In the term involving the factor XZi_iUZ,uz,_|_i , one 
just uses the bound xcoii, (...)< 1 to obtain that this term is 0{r]e), whereas in the other one, one 
may replace X"^ by At the end, 

and this last integral is 0{i]e). So 



XCoiifc+i {Xi{x, t + ri), Xi+i{x, t + r?)) dcik,k+i, 



4- Boun, 
first one: 



Bi,,,fe,i = O(e'). (5.37) 
ri'l2^=i^2,j,k,i o-nd v'iyj^i^3,j,k,i- These two cases are similar, and we treat only the 



li^U U - 1)??)) dq 



< Ce. supT]^xii-e,i]{Xi{xi, (j - l)tj)). 
Since, for any x £ ^, the number of non zero terms in the sum is 0{ne) 

0{er). 



5. Bounding 'Qi.j.k,!- 



B4,,,fc.i < % I XCoii,(x).P^'''e?7dq = 0{€^) 



(5.38) 



(5.39) 



6. Concluding the bound on fi2- The bounds ( 5.35|5T39 l show that (5.331 holds for i — 1,2,3,4 and 
1 < fc, / < TV — 1. However, as said after ( 5.27|5.3"0 l, one still needs to explain how to adapt the proof 
when one starts with one of the terms ( |5.28|5.30| l, instead of ( |5^ . If one starts from (15.281), all the 
previous arguments actually still hold. If instead one starts from (5.291, the only place where some 
non-immediate adaptation is required, concerns the bound on Bi^.^^t for t — {j — l)ri < c, where c 



is the constant introduced just after (5.341. Indeed, for x G Coll, it can be that 

{x°kixt{^,v),t),x2+i{n+i(x,ri),t)) e Ik 

for a finite and constant number of times t — {j ~ l)rj < c. Since now the factor 1/77^ in the beginning 



of right hand side of (5.341 may be replaced by a factor l/ry, one have, for such times, only the bound 

1 



Bi,j,fe,fc < 



VJq 



XCoiifc (a;) dq. = 0{e) 



But then, inserting this in (5.331, one gets ?7X], Bi,_,,fc 



0(r;e) — ©(er), since the sum contains 



only a finite and constant number of non-zero terms. The case of (5.301 is analogous. □ 
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5. The case N' < N. If iV' < iV, the arguments leading from ( |5.34[ ) to ( |5.35[ ) do not hold 
anymore, and this is the only place where a problem occurs. Indeed one assumed there that all the 
energies were different in order to make use of the diophantine condition. This problem is clearly 
superficial, since we do not care about collisions between particles having the same energy, and it can 
indeed be fixed by a slight modification of the proof. 

Let us go back to the general strategy explained in Subsection |5.1| and let us again start from a 
measure P assigning to each particle a given energy. When Duhamel's formula is used for the first 
time, in (5.2 1, instead of taking S to be the operator 'p"*-'' decoupling all the particles, one can take 



an operator which decouples only blocs of particles having a given energy. Up to an error of order e^, 
this measure is still invariant under this dynamics, and this is in fact all what we need to proceed. 
When Duhamel's formula is used for the second time, in (5.31, one still takes S to be 'p"*'"' however, 
since otherwise one could not make use of the properties obtained in the Section |6] 

It can then be checked that the proof essentially remains unchanged but that, when bounding ^2, 
one may now assume that recoUisions only occur between particles of different energies, allowing thus 
the use of the diophanitine condition. We decided however to not implement these changes explicitly 
to avoid even more heavy notations. 

6 Uncoupled dynamics 

We analyze here the dynamics when e = 0, so that all the particles evolve independently. 

6.1 = 1 particle 

Take = 1, and let ei be the energy of the only particle. Let X = (g, p) be the chain on K x P = 
M X {—61,61} solving the equations 

dq = pdt, dp = -2gdN, 

where N is a Poisson process of parameter A > 0. For a given smooth function u on K x P, let 

u{x,t) = E(« o X{x,t)), 

which, for t > 0, satisfies the equation 

dtu{q,p,t) ^ pdqu{q,p,t) + X{u{q, -p,t) - u{q,p,t)). 
It is lengthy but straightforward to check that u may be written as 

/+00 
p{q,p,q\p',t)'v{q',p')dq' 



where the probability measure p is given by 



p{q,p,q',p',t) ^* ^^'T-Pn(g,p,g',p',t) (6.1) 

^ — ^ n.' 



n< 

n = 
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with 



po{q,p,q' ,p' ,t) = X{o}ip- p')So{pt + q - q), 



Pn{q,p,q' ,p' ,t) = Xl-\p\t,\p\t]i<l' - l)X{0}{P - p') 



Pn{q,p,q' ,p' ,t) = Xl-lp\t,\p\i](l' ~ l)X{o){p + p') 



|p|t2"(f)! 



1 _ 1 ~i 
pt 



\p\t2- 



n > 2 even, 



n > 1 uneven. 



This allows us to give an expression for the kernel of the evolution operator P'^'' of the chain X 
on <1> = [0, 1] X { — ei, ei} defined by ( |3.1|3.2[ ). Given a function u on $, and given t > 0, let us write 

r°'^u{x)= ^ / f{q,p,q',p',t)u{q',p')dq'. 

The trajectories of X are obtained from these of X by reflecting them against the boundaries at q — 
and q = 1, and therefore 

f{q,p,q',p',t) = + 5Z p('?.P. 1-'?' + ^. -p'.^)- (6-2) 

fc£Z,even fcG2^, uneven 

Lemma 3. One /las this. 

1. (Doeblin's condition) There exists a > and to>Q such that f{x,x',to) > ^ for every x,x' £ <&. 

2. There exists a function g G L°°(<1?^ x R+) 5mc/i that 



f{x,x',t) = e ^'x{p(^,t)}(p') '59{i^.t)(g') + g{x,x',t). 



(6.3) 



Proof. The first part is directly established using (6.1 1 and (6.2 1. Let us show the second one. By 



(6.2 1, the first term in the right hand side of (6.3 1 directly comes from po in (6.1 1. Then, inserting 



(6.11 in (6.2l, and putting the sums over k inside the sums over n, one sees that it suffices to find a 



constant C < +oo such that, for every p G P, for every z G R and for every n > 1, one has 



E 1 + 



w 

' fcez 
|z+fe|<|p|i 



z + k 



1 - 



z + k 
pt 



if n is even, and 



L ( 

:(^n^y,(^n^y, \r>\t ^ V 



1 + 



fcez 

lz + fcl<|p|t 



z + k\ 



pt 



1 - 



z + k 
pt 



< C, 



<c. 



(6.4) 



(6.5) 



if n is uneven. Let us show (6.41 ; the proof of (6.5 I is analogous. One finds c > such that, for 
every x G [—1, 1], 

(l + a;)"/^(l^a;)"/^-i < 2(1 - x^)"/^-^ = 0(6""""'). 



But 



J_ y e-'^"(^)' = 0{n- 
\p\t ^ 



l/2x 



|z + fe|<|p|t 



On the other hand, one establishes by means of Stirling's formula that the combinatoric factor in the 
left hand side of (6.4 1 is ©(n^^^). This completes the proof. □ 
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6.2 N >1 particles 

Lemma [3] has direct implications to the case N > N' > 1. Since the particles evolve independently, 
the operator P"'* can be written as the tensor product 

V°'' =Vl...rh, (6.6) 

where, for 1 < k < N, is the evolution operator of a chain which moves only particle k and let 
fixed the other ones: 

Viu{x) ^ y2 / f{<ik,Pk,q'kiP'k,t)'u{qi,---,qk,---,<iN,pi,..-,p'k,-..,PN)dq'f,, (6.7) 

, I Jo 

with / defined by ( |6.2[ ). Now, Doeblin's condition (the first part of Lemma [3| means that the 
operators V^*'^° have a spectral gap for the || . ||i-norm. By the tensor product decomposition (6.61, 



one concludes thus that there exist constants C < +00 and c > such that, for every measure /i on 

llP^'-V-Mlle)!!, < Ce-^*||M||i, (6.8) 
with /^(l|e) defined by | |5.f 8[ ), e = (ei, . . . , ejv), and Ck the energy of particle k. 
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7 Figures 



Figure 1: The dynamical system for N — 5. Particles may collide elastically with their neighbor 
when they enter in the overlap region of size e. 



Qk+l 








Figure 2: Projection of $ on the (g^, gfc+i)— plane: the two circles represent projections of two 
possible states of the iV-particle system. Points in the triangle in the lower right corner belong to 
^ — since Qk+i < Qk — i + e there. In the region on the left of the dotted line, particle k may 
collide with particle k — 1, whereas in the region above the dotted line, particle k + 1 may collide 
with particle k + 2. 
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